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ABSTRACT 


A time series of the Loop Current in the 
Gulf of Mexico, covering an annual cycle of 
growth, spreading, and decay, has been ob- 
tained in synchronization with ERTS. Com- 
puter enhanced images, which are necessary 
to extract useful oceanic information, show 
that the current can be observed either by 
color or sea state effects associated with 
the cyclonic boundary. The color effect 
relates to the spectral variations in the 
optical properties of the water and its sus 
pended particles, and is studied by radi- 
ative transfer theory. Significant oceanic 
parameters identified are: the probability 

of forward scattering, and the ratio of 
scattering to total attenuation. Several 
spectra of upwelling diffuse light are com- 
puted as a function of the concentration of 
particles and yellow substance. These 
calculations compare favorably with experi- 
mental measurements and show that the ratio 
of channels method gives ambiguous inter- 
pretative results. These results are used 
to discuss features in images where surface 
measurements were obtained and are extended 
to tentative explanation in others. 
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This is the sixth report on project EREP 108, and it covers the month of 
December 1973. 

* . 

The first attempt at a data take was performed on 3 December. The R/V VIRGINIA 
KEY sailed the afternoon of 2 December and was on station at the intersection 
of track 62 and the Florida Keys at sunrise the next day. The line began as 
soon as the zenith angle was 80°. Continuous measurements of chloronhyll-a 
and surface radiometric temperature were taken as well as close spaced samples 
of volume scattering function and BTs. At the first (and only) station 
spectral irradiance and plankton observations were made; the cloud cover was 
751. . At the edge of the shelf the Gulf Stream peaked the seas up so that 
further stations were impossible. Further the cloud cover was increasing and 
the aircraft could obtain only two levels of the seven required. Therefore, 
at 1500Z the aircraft was radioed to abort the mission and both the ship and 
the aircraft headed back. (Later conversations with JSC revealed that the 
spacecraft would not have obtained the required data anyway) . The ship was 
readied again for a pass on 8 December but that was cancelled by JSC, 

On 11 December the P.I. gave an invited paper at the ERTS Symposium in 
Washington, D.C., entitled '"Relationships Between ERTS Radiances and Gradients 
Across Oceanic Fronts' 5 , by George A. Maul and Howard R. Gordon (co-P.I.). 

Since this lays a significant foundation for the SKYLAB mission it is ap- 
pended to this report. 

Recipients of the financial report are marked by an asterisk on the attached 
distribution list. 
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INTRODUCTION 


Remote sensing of the ocean in the visible region of the 
spectrum is being explored for a variety of reasons. Maximum 
insolation at 475 nm provides a natural energy source for 
passive sensors. A minimum in the attenuation coefficient 
for water may allow a passive measurement of some oceanic 
properties as a function of depth. Variations in the spectra 
of upwelling light at these wavelengths can be attributed to 
variations in pigment forming molecules • such as in phyto- 
plankton and to variations in the concentration of scattering 
particles such as suspended sediments. Strongly baroclinic 
currents in the subtropics and tropics can be detected by 
change in the optical properties across their boundaries when 
infrared techniques fail due to isothermal surface conditions. 
Patterns of man's activities in and on the ocean have near 
surface manifestations which change the nature of upwelling 
visible radiance. 

The social significance of oceanic observations from a 
vehicle such as ERTS has recently been summarized by Maul 
(1974a). Discussions in the present paper are limited to a 
study of the dominant feature of the Gulf of Mexico circula- 
tion, that portion of the Gulf Stream System called the Loop 
Current. This intense current transports vast amounts of 
thermal and kinetic energy into the basin through the Yucatan 
Straits. Temporal and spatial variability in the flow (figure 1) 
was studied as part of the ground truth time-series which 
followed the anti-cyclonic turning from Yucatan to the Florida 
Straits. Some dynamical results were reported by Maul (1973a) 
and will not be discussed here. As a purely descriptive 
oceanographic study, however, these data have significance in 
hurricane intensification studies such as by Leipper and 
Volgenau (1970); the investigation of the Florida east coast 
red tide Ohirphy , et al . , ’1973) used these ship tracks and a 

supporting ERTS image to document their arguments; the en- 
vironmental impact of a recent accidental jettisoning of 
cyanide canisters following a ship collision in the Gulf of 
Mexico is being studied using these data (Corwin and Richardson, 
1974). Continued reconnaissance of this current pattern is 
significant in our understanding of, and reaction to, the 
transport, which affects fishing industries, marine trans- 
portation, and public health. 


ERTS OBSERVATIONS OF OCEANIC FRONTS 

The location of the cyclonic boundary of the Loop Current 
was obtained by tracking the 22°C isotherm at 100 meters depth 
using expendable bathythermographs. The pathline of this 
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isotherm is of the order of 2 0 kilometers to the right 
(facing downstream) of the surface frontal zone (Hansen and 
Maul, 1970). This, coupled with ship observations of chloro- 
phyll-a, surface temperature, volume scattering function, and 
sea state, provides the baseline measurements to insure that 
the interpretations to follow are well founded. Figure 2 is a 
temperature-depth profile across the current along a suborbital 
track during a satellite transit; surface observations are at 
the top of the figure. This summarizes the general conditions 
of increasing temperature and salinity and decreasing chloro- 
phyll and scattering when crossing the boundary into the 
current and serves to orient them to the indicator isotherm. 
These observations confirm (Maul and Hansen, 1972) that 
changes in all these properties occur simultaneously; allowing 
an increase in confidence when a recognition decision is made 
using several variables. 

Diffuse reflectance from beneath the ocean, which is de- 
fined as the ratio of upwelling to downwelling irradiance, 
measured just above the surface, is rarely more than 0.05. 
Reflectance from the ocean's surface, which is independent 
of diffuse reflectance, can be comparable or even substantially 
larger depending on sea state, and this has been shown by Maul 
(1973b) to be a useful indicator of the current. Reflectance 
from clouds and agricultural scenes is sometimes an order of 
magnitude greater than from the ocean, even in the 500-600 nm 
wavelength (A) region. In order for the NASA Data Processing 
Facility (NDPF) to produce an image for an average scene 
radiance, the ocean signal is compressed into the lowest few 
gray scales. This is clearly illustrated in figure 3 which is 
a scanline plot across the boundary of the Loop Current from 
the multispectral scanner. The large spikes in all four 
channels (MSS 4 upper) are clouds; there seems at first glance 
to be very little change in digital number (DN) , which is pro- 
portional to radiant intensity, as a function of the sample 
number. However, careful examination shows that the average 
value of the DN at samples greater than number 950 is slightly 
larger than those before this point. It will be seen that 
this marks the transition to higher radiances caused by in- 
creased sea state in the current. 

In order to graphically display this small change over a 
two-dimensional region, computer enhancement is necessary. 
Contrast stretching for the ocean scene is accomplished by 
first studying the frequency distribution of DNs in a training 
area on the image; the area is selected to be representative 
of the ocean away from land and is large enough to be statisti- 
cally significant. As seen from figure 3, such a histogram 
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■would be strongly bimodal due to clouds. Selection of the upper 
(DN ) and lower (DN-.) .cutoff can be set at + 2 standard deviations 
about the ocean mode; in this case 9<DN<13 was used. It should 
be noted that each image will have different cutoff values and 
that the final image is quite sensitive to the range of DN chosen 
(Charnell et al . , 1973). The general expression for computing 
the stretch variable (4') is: . o 


c = 


DNu-DN 

DN U -DN| 


n 


(for DNi < DN < DN U ) 


CD 


where M is the maximum value allowed by the output device, and 
n is an arbitrary integer exponent. Equation 1 produces a. nega- 
tive of the input digital image; positive whole integers n, 
further stretch the low radiance values encountered in the ocean. 
The graphic result of using equation 1 on the data from which 
figure 3 was taken is given in figure 4. 

Figure 4 is a negative print of an area due north of the 
Yucatan Straits using MSS 5 data. Computer enhancement in this 
iamee uses only five gray scales of the 128 levels available, 
all values below DNi are set to 127 and all above DN U are set to 
0. The boundary between the resident Gulf waters (le £ t) and t 
current fright! is seen as a transition from light to dark ton-., 
respectively. Since the radiance levels m the ocean are so lov. r 
the oceanographer must resort to computer enhancement as this 

example shows. 

Figure 5 is another enhancement of the Loop Current boundary 
using MSS 5 where 7<;DN<15 and n * 2 . In this negative image of 
the western Florida Keys, water from Florida Bay extends into 
the Florida Straits and is entrained by the Loop Current. The 
current boundary in both figures 4 and 5 was_ delineated by sur 
face vessel tracks during the day of the satellite . transit . 

Notice that the current is darker in tone (higher in radiance) 
in figure 5. This is caused by the dominance of surface re 
tion due to higher sea state in the current in figure 4, as com- 
pared to higher reflectance due to particles m the Florida Bay 

water in figure 5. 

SOME THEORETICAL CONSIDERATIONS 

At this point it is useful to consider in detail the proc- 
esses contributing to the radiance spectrum N(\) at the position 
of the satellite. Solar radiation incident at the top ot tne 
atmosphere is absorbed and scattered in the atmosphere. Some 
radiation is scattered back into space without striking the 
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ocean contributing a radiance N (a) at the satellite. The rest 
(which is not absorbed) will interact with the ocean. This 
interaction can yield* upwell ing radiance above the ocean in three 
ways: 1) specular reflection from the surface; 2) diffuse re- 

flection from foam (bubbles) on or just beneath the ocean surface 
and 3) the diffuse reflection from water molecules and suspended 
particles in the water. Of these the first two phenomena are 
closely related in that they depend on the sea state, while the 
third source of radiance is essentially independent of sea state. 
The specular reflectance from the rough ocean surface can be 
computed by the methods of Cox and Munk (1954) in terms of the 
wind speed and incident radiance distribution. The diffuse re- 
flectance from white caps can be approximately accounted for by 
assuming they are "white" , Lambertian, and have an albedo of 1 
so that they contribute uniform upwelling radiance just above 
the surface given by 

H o <*> TF 

where f is the fraction of the scene covered by the white caps 
and H 0 (a) the irradiance incident on the sea surface. The third 
source of radiance, that from beneath the surface, is the most 
difficult to compute and will be discussed in detail below. We 
can write the radiance at the satellite as 

N(X)=N s (X) + y(X)N ss (X) + a(X) N d (X) ( 2 ) 

where N < ,„(A) is the contribution at the surface due to reflec- 
tion from the surface and white caps, N d is the diffuse radiance 
just above the surface due to photons that have penetrated the 
ocean, « (a) and v(a) are atmospheric transmittance factors for 
N d and N ss . a (a) and '/(a) are in general not equal since 
the radiance distribution (variation with angle) of N ss and N, 
are different. It is possible for photons to be reflected from 
the surface, backscatter from the atmosphere into the ocean, 
and scatter back into the atmosphere. Photons which do this we 
consider to be a part of N d . It should be stressed at this 
point that N d (A) in the above equation is the only source of 
radiance that contains information about conditions beneath the 
sea surface such as the concentration and composition of sus- 
pended particles and dissolved organic material. thorough 
understanding of the dependence of Nj(a) on the basic optical 
properties of the water and its constituents is required in 
order to obtain quant -it at i.v e information about the constituents 
of the ocean from measurements of N,(a). The remainder of this 
section is devoted to relating N d (AJ to these optical properties. 
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The total phase function for water and 


particle scattering 


is 
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P(.Q) - (by V P w (g) -f bppp(^)) 

(b w * b p ) 

where P w and P^ are the phase functions due to water only 
and particles only respectively. It is convenient to further 
define the forward (F) and backward (B) single scattering 
probabilities by 

B~i-F 

F=2irf”' 2 P(8) sin0 dd (6) 


and the single scattering albedo by 



Hence it is clear that 


(bpBp+b w B w ) 
( bp + b w ) 


( 8 ) 


w 0 


(bw + bp) 

(aw + op+dy + bp+bw) * 


( 9 ) 


and in general, 


0< <u 0 < | . 


It should be noted that all of the above quantities 
depend on wavelength (A) . The transfer of radiation in the 
ocean is governed by the radiative transfer equation which 
has been discussed in detail by Chandrasekhar (I960) and 
Preisendorf er (1965). Gordon and Brown (1973) have computed 
the diffuse reflectance (upwelling irradiance/incident irradi- 
ance) just above a flat homogenous ocean as a function of its 
optical properties by a Monte Carlo technique. Using a com- 
bination of the parameters which arise naturally from the 
quasi-single scattering model (Gordon 1973), Gordon, Brown 
and Jacobs (in preparation) show that the diffuse reflectance 
(Rd) can be written 

; Rd=O.I79x + 0.05I0X 2 + 0.I7I0X 3 (10) 
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In order to use equation 10, B(A) and a) 0 (A) in equations 
8 and 9 must be determined. We shall follow Tyler, Smith and 
Wilson (1972) and use for b y and B w values experimentally 
computed for sea water continuously filtered for 18 hours 
(Petzold 1972), rather than itheoretical computations of these 
quantities for pure water. We do not feel comfortable with 
this, however, on the above basis they have estimated the 
absorption spectrum of "clear natural water", and since we 
use their a w (A) spectrum this assumption is essential for 
consistency. It is thus assumed that B w = 0.1462 and 
b w = O.OOSnT 1 at 550 nm, so . 

,’b w = 0.008 [-^p] 4 m-i 

To find Brj, Petzold’ s measurements of P(0) from the Tongue of 
the 0cean P (T0T0), San Diego Harbor (SDH) ; and Off Shore Cali- 
fornia (0SC) at 530 nm have been used. Assuming the above for 
■B w and b w , the contribution from water is subtracted yielding 
the B d values listed in Table 1 (w 0 is given in the paren- 
thesis and o, 0 for filtered sea water at this wavelength is 
0.136). 


TABLE 1 

Derived Values of B p for the Indicated Locations 


T0T0 

0.0165 (.59) 
0.0158 (.25) 
0.0130 (.26) 


SDH 

0.0186 (.82) 
0.0194 (.83) 
0.0169 (.91) 


OSC 

0.00966 (.59) 
0.00836 (.55) 
0.00836 (.55) 


From Table 1 it is clear that B p not only varies considerably 
from one location to another, but also varies considerably 
for a single region. This is unfortunate in that for small 
o> , is directly proportional to <^ 0 B and w 0 and B are then 
equally important in determining R d for these cases. We 
shall assume that Table 1 gives the range. of variation of B p 
to be expected in natural waters, however this is probably 
not the case. Furthermore we will assume B p (and B^) is 
independent of wavelength. 

Since a is taken to be zero, and the a w (x) estimate of 
Tyler, Smith p and Wilson is to be used, only a y (A) remains to 
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We now turn to the problem of estimating hr)/b w or some 
related quantity from satellite observations. Since the 
above calculations are only for the case of negligible 
phytoplankton (chlorophyll) concentrations we must find 
criteria from which to choose which ERTS channels will satisfy 
this constraint. Also unless only MSS 5 is used we expect 
yellow substances to be important especially in coastal re- 
gions where river runoff, etc., may be considerable. Further- 
more there is the additional problem that the reflectance 
depends on B p which is also unknown, as -will as bp/b w . We 
have developed a method which partially overcomes some of 
these problems. Considering the radiances observed in MSS 4 
and 5, we have approximately 


N 4 =a 4Nd 4 + r4Nss 4 +N S4 

N 5 =a 5 N d 5 + r5Nss 5 + Ns 5 C11) 

where and 7 i (i = 4, 5) are the fractions of Ndi and N ss ^ 


(measured at the sea surface) that reach the sensor. It is 
assumed that and are constant over an ERTS Frame (if 
they vary in a known way, their influence is easily accounted 
for and will not be discussed further). Now from the theory 


Nd 4 = k 4 




( 12 ) 


where and k^ are essentially independent of bp/b w , but 
depend directly on Bp. Taking the horizontal gradient, 
of equation 11 using equation 12 we find 

b, 


7 H Ni = k |ai V H ^j + yi V H N ss 

w 


(13) 


since N S | is nearly constant over a frame (VN S . = 0). The last 
term in equation 1 13, V^N-s^, is the horizontal gradient of the 
reflected radiance from the sea surface. This is nearly zero 
everywhere except where the sea state changes dramatically 
with horizontal distance (for example in figure 4) . Nearly 
everywhere on the frame then, V^N SS ^ = 0, so 

7 H N 4 _ k 4 a 4 

V H N 5 k 5 a 5 
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properties of fhe individual components with wavelength is 
known . 

As an example of this, the spectra given in figure 12 
were integrated over the MSS 4, MSS 5, and MSS 6 filter 
functions. A series of numerical tests were then made of 
ratios, differences, and sums to see if the three water types 
could be distinguished at the sea surface. It was quite easy 
to distinguish on the basis of such calculations between the 
Gulf Stream waters and the coastal waters, and between the 
coastal waters and the plankton bloom; however it was not 
possible to distinguish between the Gulf Stream and the 
plankton bloom. This suggests, as the theory implies, that 
for these data the ratio test (MSS 4/MSS 5) is not likely to 
be successful in specifying the chlorophyll -a concentration. 
Further, since the sea surface component, N^ s , (a) 5 spectrally 
alters N (A), numerical tests (ratios, differences etc.') are 
invalid indicators of oceanic properties. 

Probably the most efficient method of determining the 
concentration of the constitutents in the ocean will be to 
compare theoretical and experimental spectra, adjusting the 
constitutent concentrations in the theoretical spectra until 
agreement is found. This of course requires a basic under- 
standing of the optical properties of the constitutents which 
can be derived only from careful in ■situ and laboratory experi- 
ments. It seems that at the present time much energy and 
money is being expended to try and use optical methods to 
locate and study materials with nearly unknown optical pro- 
perties suspended or dissolved in a medium with only poorly 
known optical properties. This must be overcome before 
significant progress can be made. 

A final interpretation of oceanic observations from ERTS 
is given in figure 13 which is a computer enhanced MSS 6 
negative image of southeastern Florida. The dark lineation 
parallelling the coast in the upper portion of the image is a 
zone of high reflection caused by locally increased N ss along 
the edge of the Florida Current. The increase in surface 
reflectance is probably caused by surface wave interaction with 
the cyclonic boundary and is not bottom influence. This further 
explains the edge effect and supports the discussion on the 
local dependence of v^N on v h^ss only given above. Another 
example of the dominance of N ss is the bright slick areas 
(low N) off the Virginia Key sewer treatment plant. This is 
probably caused by the dampening of the glitter causing 
capillary waves in' the oil film associated with the organic 
slick. The slick, which has drifted south past the popular 
Key Biscayne beaches, offers an explanation of the narrow 
lineation off shore in the Florida Straits; a passing oil 
tanker heading south which is pumping her bilges would cause 
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where the least squares values for g=0.38 and h=0.00, and the 
linear correlation coefficient, r , is 0.66. The coefficient of 
determination, r 2 , which is the ratio of the explained varia- 
tion from the mean (by the least squares line) to the total 
variation, is 44%. Physically this means that this image may 
not be useful for determining particle concentrations because 
B d , w 0 (i.e. a v ) t or b p /b w may be changing. The extent to_which 


this holds true in natural waters is unknown; this will be 
extensively tested in the New York Bight area where turbidity 
measurements are being made by NOAA vessels concurrent with 
ERTS transits. 


The data of figures 7 and 8 can be compared to measure- 
ments of upwelling irradiance. Data given in figure 12 were 
observed using a %-meter Ebert spectroradiometer from 3 meters 
above the surface; these observations were made during the 
time frame of figure 5 and represent the water types shown in 
that image. All spectra were carefully selected to represent 
the same downwelling irradiance, sea state, sun angle, cloud 
cover, and absence of bottom influence. Specular reflection 
due to waves were minimized by preselecting ten spectra with 
similar shapes. After digitizing the records, averages and 
standard deviations (<0 were computed at each wavelength; if 
values exceeded the average by 1 they were rejected and a 
new mean computed. Absolute values of the spectra are traced 
to NBS through the 2 -meter integrating sphere at NASA’s 
Goddard Space Flight Center. 


As the chlorophyll-a concentration increases, the spectral 
peak shifts to longer wavelengths in similarity with the com- 
puted data in figure 8. This is not to say that the chloro- 
phyll-a and the yellow substance produce quantitatively similar 
results; in general they will not. What it does imply is that 
the broad absorption bands in pigments such as chlorophyll and 
yellow substance at shorter wavelengths (<600 nm) produce 
qualitatively similar effects in the spectra. When other fac- 
tors are equal, increasing the concentration of pigments will 
cause a decrease in the radiance in MSS 4 only . In nature 
however, such as in a plankton bloom or in a river plume, 
increased amounts of chlorophyll are normally accompanied by 
increased particle concentrations (organisms which contain 
the chlorophyll) as well as yellow substances (decay products 
of the plankton) , and increased amounts of ’’white” ^ suspended 
particles in the river plume case. The radiances in MSS 4 
and 5 will in these cases vary in a maimer whieh frill depend 
again on the relative concentration of the constitutents , and 
since a and b vary almost independently, the reflectance 
signature, which depends on w o(A), and B(A) is not unique and 
can be unraveled only when the variation of the optical 
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a similar feature on the image. Thus ERTS could be useful in 
patrolling coastal waters for such illegal acts which affect 
the nearshore water quality. 


CONCLUSIONS 

spectral properties of the oceanic front associated 
with the Loop Current have been studied by ship and satellite 
observations and by radiative transfer theory. It is seen 
that computer enhancement is required to extract useful 
information from the ERTS data for the ocean scene. The 
current boundary can be detected by changes in the surface 
reflectance, N ss , as well as the diffuse reflectance, Nar from 
below the surface; N^ however is dependent on both E and 
and thus the spectral interpretation of ocean color requires 
surface truth measurements for meaningful results. Particle 

whi ? h can delineate currents, can be estimated 
m MSS 5 if the ratio t H N 4 /t„N 5 is reasonably constant over a 
scene. MSS 4 is. strongly influenced by yellow substance, and 
particle estimation based on these data are invalid in many 
z °hes. Water mass identification using ratios or 
differences of MSS 4, 5, and 6 data have no validity in either 
theory or ^ observation in the Gulf of Mexico. Finally, it must 
be emphasized that the sp.e.ctrum o.f upwe.lling radiance just 
above the surface is a function of both Nd and and that 
N ss frequently dominates. 
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Figure 1. Time series of the Loop Current in the Gulf of Mexico, August 1972- 
September 1973. The pathlines are the location of the 22°C isotherm 
at 100 meters depth. Each cruise was synchronized with ERTS passes 
over the area every 36 days. The indicator isotherm was located by 
expendable bathythermographs from a surface vessel; this isotherm is 
approximately 20 kiiometers to the right of the cyclonic edge facing 
downstream. 
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Figure 2. Cross section of the Loop Current taken by ship along a suborbital 
track on the day of ERTS transit. Section shows typical increase 
in salinity and temperature (dashed-radiometric, solid-bucket) when 
crossing into the current. Similarly chlorophyll-a and volume 
scattering function at 45° (436nm light) decrease. The temperature- 
depth section shows the relation between the 22°C at 100 meters and 
the surface frontal zone (at ca. 1330 GMT). Horizontal scale is 5 
hours travel time at 9.8 knots or 90 kilometers. 
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SCAN SPOT NUMBER 

Figure 3. ERTS scanline plot across the Loop Current front. Top scanline is MSS 4, next MSS 5, MSS 6, and 
MSS 7 on the bottom. The large energy spikes are clouds. At scan spot number 950 there is an 
increase in the average value of the digital number of 1 or 2 ; this marks the cyclonic edge of 
the current. 
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Figure 4. Negative print of computer enhanced (9<DN<13; n - 1) MSS 5 image of 
the cyclonic boundary of the Gulf Loop Current. Surface vessel track 
confirmed the location of the current to be the darker shade (higher 
radiance) region on the right hand side of the image (ERTS ID 1065- 
15411). Scanline plot in Fig. 4 horizontally passes through the 
middle of the scene. Horizontal distance across the image is 90 
kilometers. 
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Figure 5, Negative print of computer enhanced (7<_DN<_1 5; n = 2) MSS 5 image of Marquesa Key and Key West (ERTS 
ID 1099-15293). Change in radiance southwest of Marquesa from' dark to light marks the ship-located 
boundary between the higher intensity Florida Bay water and the lower intensity Gulf Stream. Bottom 
depth is in excess of 100 meters and thus does not contribute to the radiance. Horizontal distance 
across the image is 135 kilometers. 
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Figure 6. Computed reflectance in percent at the sea surface as a function of the ratio of the particle scat- 
tering coefficient to the water (only) scattering coefficient for 550nm, The value, of the fraction 
of backscattered light due to particles (Bp) for each curve is given on the right hand side. 
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Figure 7, Computed reflectance in percent at the sea surface as a function of 
wavelength for various values of the ratio of the partiple scattering 
coefficient to the water scattering coefficient. Note the wavelength 
dependence in these spectra of changing the particle concentration 
only. 
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Figure 8. Computed reflectance spectra, in percent, at the sea surface. The 
ratio of the particle to water scattering coefficients ( b p /b w ) and 
the fraction of backscattered light are 128 and 0,0165 respectively. 
Values of the absorption coefficient due to yellow substance (a v ) are 
listed on the left hand side. Note the shift of the spectral peak to 
longer wavelengths with increasing a . 
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Figure 10. Computed reflectance at the sea surface, in percent, as a function of the particle to water scat- 
tering coefficient ratio, integrated over the response function of the MSS 5 filter (0.6-0,7ym). 
Valuesof the absorption coefficient due to yellow substance (a ) are listed along the right hand 
side. Note the linearity is not strongly dependent on a y but only on the fraction of backscattered 
light, B_ . 




Figure 11. Contrast stretched (4<DN<1 2; n = 2) negative image of the ocean area offshore of Cape Hatteras 
(ERTS ID 1132-15042). TFe Gulf Stream can be seen as the bright area to the south of the en- 
trained sediment from the coastal estuaries. The least squares fit of eq. (17) was done along 
a scanline north of the Cape and extending from nearshore, through the suspended sediment and 
into the current. Extensions of this plume were observed for 150 kilometers further east on 
other ERTS images. Horizontal distance across the image is 135 kilometers. 
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Figure 12, Observed upwelling spectral irradiance in the Gulf of Mexico, 

November 1972, The three spectra represent typical observations 
during the time series, and show the shift of the dominant wave- 
length to larger values with increased surface chlorophyll-a. 

The volume scattering function, 645, is for blue (436nm) light. 
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Fiqure 13, Negative print of computer contrast stretched (7<DN<15; n = 2) F.RTS 
image of South Florida (ERTS ID 1026-15230). The Florida Current 
can be seen as a line of dark lineation parallel to the coast; bot- 
tom is essentially invisible in this MSS 6 scene. A ship can be 
seen by its characteristic V-shaped wake just offshore of Miami 
Beach. Possibly the Virginia Key sewer outfall area ran be observed 
by its low reflectance due to an organic slick. Horizontal distance 

across the image is 90 kilometers. 
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